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The Higgs boson production pp (pp) —> bh + X via bg bh at hadron colhders, which may be an 
important channel for testing the bottom quark Yukawa coupling, is subject to large supersymmetric 
quantum corrections. In this work the one-loop SUSY-QCD corrections to this process are evaluated 
and are found to be quite sizable in some parameter space. We also study the behavior of the 
' corrections in the limit of heavy SUSY masses and find the remnant effects of SUSY-QCD. These 

' remnant effects, which are left over in the Higgs sector by the heavy sparticles, are found to be 

. so sizable (for a light CP-odd Higgs and large tan /3) that they might be observable in the future 

' experiment. The exploration of such remnant effects is important for probing SUSY, especially in 

case that the sparticles are too heavy (above TeV) to be directly discovered in future experiments. 
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' I. INTRODUCTION 

> 

Searching for the Higgs boson is the most important task for the Fermilab Tevatron coUider and the CERN Large 
■ Hadron CoUider (LHC). Among various Higgs production mechanisms, those induced by the bottom quark Yukawa 
coupling are particularly important because in some extensions of the Standard Model (SM) such a coupling could 
be considerably enhanced and thus the production rates can be much larger than the SM predictions. The Minimal 
^ Supersymmetric Standard Model (MSSM) [1] serve as a good example of such extensions, where the coupling of 
. the lightest CP-even Higgs boson (denoted by h) to the bottom quark is proportional to tan/3 [2] and thus can be 
' significantly enhanced by large tan/3. 
Oh, In the production channels of the Higgs boson via its coupling to the bottom quark , the process pp (pp) bh + X 
via bg — > bh was recently emphasized in Ref. [3] . The advantage of this process over the production via bb ^ h [4] , the 
^ , dominant production channel via the bottom quark Yukawa coupling, is that the final bottom quark can be used to 
' reduce backgrounds and to identify the Higgs boson production mechanism [5]. And compared with the production 
via gg,qq — > hbb [6], the production rate of pp (pp) — ^ bh + X is one order of magnitude larger. So the production 
pp (pp) bh + X may be a crucial channel for testing the bottom quark Yukawa coupling. 

If the MSSM is indeed chosen by Nature, then the prediction of the cross section for the production (pp) bh+X 
[7] must be renewed with the inclusion of SUSY quantum corrections because, like the process of the charged Higgs 
boson production pp (pp) — tH~ + X [8,9] and the relevant Higgs decays [10-12], the SUSY quantum corrections to 
this process may be quite large. In this work we study the one-loop SUSY-QCD corrections to this process, which is 
believed to be the dominant part in the SUSY corrections. 

It is well known that the low-energy observables in the MSSM will recover their corresponding SM predictions when 
all sparticles as well as Ma (the mass of the CP-odd neutral Higgs boson) take their heavy limits. If only some of the 
masses take their heavy limits, e.g., all sparticles are heavy but Ma is light, then large remnant effects of SUSY may be 
left over in the physical observables of the Higgs sector. The study of these remnant effects will serve as an important 
probe for those heavy SUSY particles [11]. Such kind of study will be performed for the process pp (pp) ^ bh + X in 
this work. After deriving the SUSY-QCD corrections to this process, we will examine the behavior of the corrections 
in the limit of heavy SUSY masses. When Ma is light and the sparticles take their heavy limits, we find the large 
remnant effects left over by the SUSY-QCD in such a Higgs production process. 

This paper is organized as follows. In Section II we present our strategy for the calculation of the one-loop SUSY- 
QCD corrections to the process pp (pp) bh + X . In Section HI, we scan the parameter space of the MSSM to 
estimate the size of the SUSY-QCD corrections. In Section IV, we study the behaviors of these corrections in the limit 
of heavy SUSY masses. The conclusion is given in Section V and the detailed formula obtained in our calculations 
are presented in the Appendix. 
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II. CALCULATIONS 



At high energy hadron colhders, the incoming 6-quark is generated from gluons spHtting into nearly colUnear bb 
pairs. When one member of the pair initiates a hard-scattering subprocess, its partner tends to remain at low pT and to 
become part of the beam remnant. Hence the final state typically has no high-pT 6-quarks. When the scale of the hard 
scattering is large compared with the 6-quark mass, the 6-quark is regarded as part of the proton sea [13]. However, 
unlike the light quark sea, the 6-quark sea is perturbatively calculable. If the scale of the hard scattering is /i (for the 
scale we use fl to distinguish from the SUSY parameter /i), the 6-quark distribution function h{x,jl) is intrinsically 
of order Q!s(/i) log(/i/m(,). As jl approaches rrih from above, b{x,jj,) 0; while as p, becoming asymptotically large, 
as{jl)log{jl/mb) approaches order of unity and one needs to sum terms of order a"{fl) log"(/x/m6) into the 6-quark 
distribution function to yield a well-behaved perturbation expansion in terms of [13]. In this case, the 6-quark 
distribution function becomes of the same order as the light partons. The main uncertainty of 6-quark distribution 
function comes from that of gluon distribution function which is about 10% [14]. 

The subprocess gb bh occurs through both s-channel and t-channel shown in Fig. l(a, 6). The spin- and color- 
averaged differential cross section at tree-level is given by 

da" _ asifi) f grrib{fi}\^ f sma\^ 1 rri j-^ + [.^ + i - j u'if 



di 24 \ 2mw ) \ cos (i J si 

where s and t are the usual Mandelstam variables, as(At) is the running strong coupling, and nibijl) is the running 

bottom quark mass [4]. a represents the mixing angel between the two CP-even Higgs boson eigenstates and is 
defined by tan/3 = V2/V1 with vi^2 denoting the vacuum expectation values of the two Higgs doublets [2]. In Eq.(2.1), 
we use the MS running mass of the 6-quark rather than the pole mass to take into account large QCD logarithm 
corrections to the vertex hbb [15]. The SM prediction of the cross section is recovered when setting | sina/cos/3| = 1 
[3]. Throughout the calculations we neglect the 6-quark mass except in the 6-quark Yukawa couplings. 

The one- loop Feynman diagrams of SUSY-QCD corrections are shown in Fig. l(c-r). In our calculations wc use 
dimensional regularization to control the ultraviolet divergences in the virtual loop corrections. For the renormalization 
of strong coupling constant Qs, we employ the MS scheme [16]. As to the hbb Yukawa coupling, at one loop level to 
0{as) it is given by 

9hbb = ghbb + Sg'iS,'' + ^gl^b"""^ (2-2) 

where ghbb denotes the one-loop coupling, g^^b is tree-level coupling, Sg'^^^ is the radiative correction from pure 
QCD [15], and ^^f^*"^ is the one-loop SUSY-QCD contribution [11]. In determining Sg^^^ + ^g^^b^^ , one needs 

the counter- term of the vertex hbb, whose general from is given by ghbbi 1 — ) with 6mb being the 

counter-term of the 6-quark mass defined by m[J = nib + Snib {m!^ is the bare mass) . Srub is determined by requiring 
rub to be the pole of the one-loop corrected 6-quark propagator [15,11,17]. One major difference between Sg'j^^i^'^ and 
Sg^bb^^ is that the former contains large logarithms as log of 0(1) and thus one needs to introduce MS running 
mass rubip) to absorb leading logarithms Q!"log(^)" [15]. An extensive discussion about this issue in the MSSM 
was provided in Ref. [17]. 

The one-loop SUSY-QCD contribution to the amplitude of gb — * bh can be written as 

5M = !^-^f^w(p2)(Ci7''PL + C27''Pfl + /kPL + ^47" APr + CnptPL + C^p'^Pr 
loTT"' 2mw cos fi 

+C7p^ /cPl + C^pt /kPR + CsP^Pl + CioP^Pfl + Cnp^ /:Pl + Ci2P^ /:Pfl)u(pi)e^(fc), (2.3) 

where Pl,r = (1 T 75)/2, = A"/2 with A° being the Gell-Mann matrices, and fc, pi and p2 are the momentum of 
the incoming gluon, incoming 6-quark and the outgoing 6-quark, respectively, gs and rub should be understood as the 
nmning ones in Eq.(2.1). The coefficients Ci arise from the loops and are given explicitly in Appendix A. We have 
checked that all the ultraviolet divergences canceled as a result of renormalizability of the MSSM. 
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FIG. 1. Feynman diagrams of gb bh with one-loop SUSY-QCD corrections: {a,b) are tree level diagrams; (c — e) are 
one- loop vertex diagrams for s-channel; {f ~k) are one-loop vertex diagrams for t-channel ; (/ — n) are the box diagrams; (o — r) 
are self-energy diagrams. 



The differential cross section of gb — > bh with one-loop SUSY-QCD corrections is given by 

da _ da° d(Ao-) 



dt 



dt 



dt 



(2.4) 



where the first term is the tree-level result given in Eq.(2.1) and the second term is the one-loop SUSY-QCD corrections 
given by 



d(Ag) 
di 



sma 



2mw 



1 



48 An \cos/3 



+ (C9 + Cio)- 



The cross section of gb bh is then given by 



dt 



da 



(2.5) 



(2.6) 



where i^^ax = and = ' -^^ 



s + mj^ with m/j denoting the Higgs mass. In order to avoid coUinear divergence in 



Eq.(2.6) and to enable the outgoing 6-jet to be tagged by silicon vertex detector at the Tevatron and the LHC , we 
require the transverse momentum of the outgoing 6-jet to be larger than 15 GeV and apply a rapidity cut |?7b| < 2.5 
for the LHC and |?7f,| < 2.0 for the Tevatron. 

The total hadronic cross section for pp [pp) ^ bh + X can be obtained by folding the subprocess cross section a 
with the parton luminosity 



cr(s) ~ [ dr —— a{s 



st), 



(2.7) 
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where tq = m^/s and s denotes the pp (pp) squared center-of-mass energy. dL/dr is the parton luminosity given by 

^=J -[f^ix,Q)f!{r/x,Q) + {g^b)], (2.8) 

where fj^ and arc the 6-quark and ghion distribution functions in a proton, respectively. In our numerical calcula- 
tions, we used the CTEQ5L parton distribution functions [18]. We did not distinguish the factorization scale Q and 
the renormalization scale jl, and assumed jl = Q = ruh- The scale dependence of our results will be briefly discussed 
in the proceeding section. 

The process pp {pp) hh + X has been extensively studied [3] in the framework of the Standard Model. Its 
cross section is found to be at the order of 1 fb for the Tevatron and 100 fb for the LHC, and the next-leading-order 
(NLO) QCD correction can enhance the production rate by 50% ~ 60% for the Tevatron and 20% ~ 40% for the 
LHC [3] , depending on the applied cuts and the Higgs boson mass. We will incorporate such QCD corrections in our 
calculations for the production rate a/asM- 



III. NUMERICAL RESULTS 



In this section we will perform a scan over the SUSY parameter space to show the possible size of the SUSY-QCD 
corrections. Before performing numerical calculations, wc take a look at the relevant parameters involved. For the SM 
parameters, we took mw = 80.448 GeV, mz = 91.187 GeV, to* = 174.3 GeV, rhbirhb) = 4.2 GeV [19], sin^ 0^ = 0.223 
and as{mz) = 0.118 . We used the one-loop QCD running as{jl) and mb{fi). 

For the SUSY parameters, apart from gluino mass, the mass parameters of sbottoms are involved. The sbottom 
squared-mass matrix is [2] 

„ / m? mtXh \ 



where 



m, 

h 



mf^ = -Km^ -Hm|(4 - Qb sin^ e^vr) cos(2/3), (3.2) 
^ = m|, + ml + rrizQb sin^ Ow cos(2/3), (3.3) 
Xb = Ab- utan p. (3.4) 



Here and are soft-breaking mass terms for left-handed squark doublet Q and right-handed down squark D, 

respectively. Ab is the coefficient of the trilinear term HiQD in soft-breaking terms and ^ the bilinear coupling of 
the two Higgs doublet in the superpotential. /| = —1/2 and Qb = —1/3 are the isospin and electric charge of the 
b-quark, respectively. This mass square matrix can be diagonalized by a unitary rotation 



6l \ ^ / cos Ob - sin Ob\ / h 
In J \ sin 6b cos 6b / \ 62 

and consequently db and the masses of physical sbottoms &i,2 can be expressed as 



(3.5) 



tan 206 = ^"^^^^ (3.6) 

m?^ = m?^ cos^ ^6 -|- 2mbXb cos Ob sin Ob + m?^ sin^ Ob, (3.7) 
m?^ = sin^ 6b — 2mbXb cos Ob sin 6b + m?^ cos^ 6b- (3.8) 

Prom Eqs.(2.1,2.5,2.6) we know that the cross section also depends on the Higgs mass, a and /3, which can be 
determined at tree level by tan /3 and the CP-odd Higgs mass [2] . Noticing the fact that both the mass and the 
mixing angle receive large radiative corrections when the SUSY scale is high above rrit [20] , we used the loop-corrected 
relations of Higgs masses and mixing angle [21,22] in the computation of cross section. In our calculation, we used 
the program SUBHP0LE2 [21], where two-loop leading-log effects of the MSSM are incorporated in the Higgs masses 
and the mixing angel, to generate nih and a needed for our computation. The input parameters for this program are 
the mass parameters in the top sqaurk and sbottom sector, and M^, tan/3 and the heavier chargino mass m^^. 
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We found that the usage of the loop-corrected relations of Higgs masses and mixing angle is indeed necessary. 
Comparing with the results obtained by using tree-level relations for Higgs masses and mixing angel, the size of 
SUSY-QCD corrections by using the loop-corrected relations is generally magnified from 30% to 200%^. We also 
checked that this conclusion is also valid for the SUSY-QCD correction to the Higgs partial width r{h bb). 

Now we know the relevant parameters are 



mQ,m^,m^,At,b, mg,m^, fi, Ma, tan/3, 



(3.9) 



where is the soft-breaking mass term for right-handed top-squark and At the coefBcient of the soft-breaking 

trilinear term H2QU. To show the main features of SUSY effects in pp (pp) —^bh + X, we performed a scan over this 
ten-dimensional parameter space. In our scan we make no assumptions about the relations among these parameters 
to keep our result model-independent, but restrict the parameters with mass dimension to be less than 2 TeV. In 
addition, we consider the following experimental constraints: 

(1) /i > and tan/3 in the range 5 < tan /3 < 50, which seems to be favored by the muon g — 2 measurement [23]. 

(2) The LEP and CDF lower mass bounds on Higgs, gluino, stop, sbottom and chargino [24,25] 



ruh > 114 GeV,mi^ > 86.4 GeV, m^^ > 75.0 GeV, rrig > 190 GeV, mf'* > 67.7 GeV, 



(3.10) 



where rn^^^ is the mass of the lighter chargino. 
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2. The scatter plot of a''^ versus Ma- 



It would be interesting to first scan over the allowed parameter space to figure out how large the production rate is 
enhanced in the MSSM. In Fig. 2 we present the tree-level cross section relative to the SM prediction with the same 
Higgs mass. This ratio is independent of collider energy, but is dependent on the SUSY mass parameters since we use 
the loop-corrected relations of the Higgs masses and the mixing angle (see the second paragraph of Sect. II) . From 
Fig. 2 one sees that the production rate in the MSSM can be significant larger than the SM prediction for a light Ma] 
while for a heavy Ma the MSSM prediction approaches to the SM value. This character was first noticed in [26] and. 



^The main reason for such an enhancement is that for a large SUSY scale, the dominant term of the SUSY-QCD correction 
is proportional to cot a -f tan /3 (see for example, Eqs.(4.2, 4.4, 4.6, 4.9) in the proceeding section), whose value obtained by 
using the loop-corrected relations of the Higgs masses and the mixing angle is generally larger than that by using the tree-level 
relations. 
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as a result of this character, distinguishing the Hghtest MSSM Higgs boson from the SM Higgs boson in the large 
Ma limit will be very difficult. When SUSY-QCD corrections are added, this character remains unchanged for heavy 
sbottoms(see following discussions). From Fig. 2 one also finds that there exists the possibility (although very rare) 
that the MSSM cross section is suppressed to be below the SM value [27]. In this case, the SUSY-QCD corrections 
will play a more important role in Higgs phenomenology at colliders [28] . 

Now let us scan over the allowed parameter space to show the possible size of the SUSY-QCD corrections relative 
to the tree-level value. In our numerical evaluation, we found the relative correction is insensitive to collider energy. 
The difference of the results between the LHC and the Tevatron is at the level of parts per mill. In Fig. 3 we present 
the SUSY-QCD corrections to the cross section. 
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FIG. 3. The scatter plot of the SUSY-QCD correction IS.a/cr''^ versus Ma for the Tevatron. The difference of the results 
between the LHC and the Tevatron is at the level of parts per mill. 



Fig. 3 manifests three features of SUSY-QCD corrections for large tan (3. The first one is that the correction size is 
enhanced by tan (3 and thus can be quite large. The second one is that for Ma lighter than 500 GeV, the correction 
tends to be negative. The third one is that for large Ma, the correction may be positive and the maximum value 
seems to be independent of the value of Ma- These features can be explained as follows. 

For the correction size larger than 2%, the dominant contribution of the SUSY-QCD correction is from the loop 
corrections to the vertex hhh and the contribution of the box diagrams is much smaller for the parameters satisfying 
the constraints in the paragraph following Eq.(3.9) ^. As a result, the correction behaves like (which is similar to the 
SUSY-QCD correction to the vertex hhh in case of heavy sbottoms [11]) 







(3.11) 



where Mew and M-j^ denote the electroweak scale and the typical mass of sbottoms, respectively. Ci and C2 are 

^ and Ah, but independent of Ma- It is found that in general Ci is negative and C2 is 



functions of mr , 

61 ' 



positive and either Ci or both Ci and C2 are enhanced by large tan/3. For a light Ma compared with to^, the first 
term of the RHS in Eq.(3.11) is dominant and hence the cross section tends to be negative. While for a large Ma , 
the second term is dominant and the cross section tends to be positive. So the behavior of Eq.(3.11) can explain the 
features of Fig. 3. 

From Fig. 3 we noticed that in some corners of parameter space the one-loop SUSY-QCD contributions to the cross 
section are comparable or even larger than the tree-level result and consequently, one must consider higher order 
corrections. In such cases, it is important to sum over the terms (^"( mJusy ^'^ orders of perturbation theory by 
using an effective Lagrangian approach [17,29]. 



In the large limit of SUSY mass parameters discussed in the proceeding section, we have checked that, even for the correction 
size far smaller than 1%, the dominant contribution still comes from the corrections to the vertex hbb. 
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FIG. 4. The SUSY-QCD correction Sa/a° and the cross section a/asM versus the mass of the produced Higgs boson. Solid 
curves are for the LHC and the dashed for the Tevatron. For Aa/a" , each sohd curve overlaps with the corresponding dashed 
one due to the tiny difference. 

Next we study the dependence of the SUSY-QCD correction 6a ja^ and the cross section normahzed by the SM 
prediction, (j/<tsm, on the mass of the produced Higgs boson, which can be directly compared with experiment 
results and hence is much informative. In such a study we assumed a common value {Msusy) for all input SUSY 
mass parameters and, considering the fact that is insensitive to Ma for Ma > ISOGeU [22], we fixed the value of 
Ma- Then through varying the value of Msusy, we obtain the different mass value of the Higgs boson. 

The dependence on mh is illustrated in Fig. 4 for tan /? = 50. (Note that ruh can vary only in a small range since 
it is stringently upper bounded in the MSSM.) In this figure and also in the following figures showing a/asM, we 
also incorporated the conventional QCD corrections [3] into both the MSSM and the SM cross sections. As pointed 
out earlier, the SUSY-QCD correction 5a /a^ is not sensitive to collider energy. The difference of the results between 
the LHC and the Tevatron is too small to be visible, as shown in the upper part of Fig. 4. But for a/asM the 
difference between the LHC and the Tevatron is visible since the QCD corrections are significantly different for these 
two colliders [3]. 

We also studied the dependence of the production rate on the renormalization scale fi. (As pointed out earlier, we 
assume that the factorization scale is equal to the renormalization scale.) We found that such a scale dependence is 
significant in some parameter space. For example, for the LHC with Ma = 300 GeV and rrih — 120 GeV, the ratio 
a{fl)/a{mh) is 0.93 for p, = mh/2 and 1.03 for jl = 2mh. Such an uncertainty is comparable with the uncertainty 
from the 6-quark Yukawa coupling {fhf,{fhh) = 4.2 ± 0.2) and the partron distribution fimction (about 10%). 



IV. BEHAVIOURS OF SUSY-QCD CORRECTIONS IN DECOUPLING LIMITS 

To study the behaviors of the SUSY-QCD correction in the large limit of SUSY mass parameters, we consider 
four typical cases as in Ref. [11] where the decoupling property of SUSY-QCD correction to the coupling of hbb is 
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analyzed. To qualitatively understand the feature of each case, we present the approximate formula in the limits, but 
in practical numerical calculations we use the complete one-loop expressions. 

About the inputs of the SUSY parameters, there are several differences between our work and Ref. [11]. The first 
one is that in [11] the tree- level relations for the Higgs masses and the mixing angle were used, but in our calculations 
we use the loop-corrected relations. As discussed earlier, using the loop-corrected relations leads to a significantly 
different correction. The second one is that in our analysis we considered the experimental bounds in Eq.(3.10). This 
will rule out some parameter space which have been considered in Ref. [11]. The third one is that in cases B and D 
we also require At^b to be large since large At^t is favored by the Higgs mass bound. 
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FIG. 5. The SUSY-QCD correction Aa/a° versus Ms in Case A. For each fixed value of Ma, the solid curve (for the LHC) 
overlaps with the corresponding dashed one (for the Tevatron) due to the tiny difference. 

(1) Case A: All SUSY mass parameters except Ma are of the same size (collectively denoted by Ms) and tend to 
heavy, i.e.. 



In this case the SUSY-QCD correction behaves like 



Ah - At 



fx - Ms. 



(4.1) 



Act 



2a, 
Stt 



-(tan/3 + cot a) — cota( 



12Af| 



ml tan^ ^ tan f3 cos f3 sm{a + (3) m 



2M| 



Ml 



(4.2) 



where the first term in the RHS corresponds to the first term in Eq.(3.11)^ and the rest corresponds to the second 
term in Eq.(3.11). The striking feature of this case is that for very large Ms, the correction approaches a nonzero 



In the MSSM the tree-level relation for Higgs masses and mixing angle predicts the following relation: cot a 
jY- tan P cos 2/3 -I- O(-r^) and at loop level, the gap between cot a and — tan /3 is generally enlarged. 



- tan (5 
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constant, and this remnant effect of SUSY-QCD corrections is enhanced by large tan/?. This feature is illustrated in 
Fig. 5. From Fig. 5 one also finds that the SUSY-QCD correction in this case is negative and sizable * for a light Ma 
and a large tan /3. 

In Fig. 6 we show the loop corrected cross section normalized by the SM prediction. From this figure we see that for 
Ma of several hundred GeV, although the tree-level cross section in the MSSM can be reduced by large SUSY-QCD 
corrections, an enhancement of several times over the SM prediction can still be expected due to the fact that the 
tree-level 6-quark Yukawa coupling in the MSSM is significantly enhanced by large tan/3 for light AIa- This large 
enhancement shows a very weak dependence on Ms- So we can conclude that up to the next-leading order, a light 
Ma is still able to make the MSSM cross section larger than the SM prediction. 
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FIG. 6. a jasM versus Ma in Case A for the LHC (solid) and the Tevatron (dashed). 

(2) Case B: Mq M^y, Mj^ and At^ (collectively denoted as Ms) is much larger than /i, rrig and Ma (collectively 
denoted as M ), i.e., 



Mr 



,U,D 



At,b ^ Ms ^ rrig fi r-. Ma ^ M. 



(4.3) 



In this case the SUSY-QCD correction behaves as 



2as 
Stt 



2M2 mlM^ ,Ms 

(tan 13 + cot a) - . (— -Hcota) 



6M| ' M 



m| cos f3 sin (a + f3) Ms ^ «^ SAf^ ' 

% (1 — ( tana) ^) 

2A/| sin a ^ ^ M ^' Mj ' 



(4.4) 



*If the correction is too sizable (say exceed 50%) , higher order corrections are also important and need a proper treatment 
[17,29]. 
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From this expression we see that in the large Ms limit, the SUSY-QCD corrections decouple rapidly as M'^/Mg and 
the decoupling behavior is slowed down by large tan/3. The characters of this case are shown in Fig. 7. So we see that 
even with a fixed light Ma, the process still does not have remnant SUSY-QCD effects if the gluino mass and are 
also kept light. 



-0.03 



b -0.06 



b 



b 



b 



-0.09 



-0.12 



3.5 



2.5 



4IW 



400 




800 



800 



Case B 



M,= M~ = [i =200 GeV 



'00 



SO 



tan fi= 30 
fan P = 8 



1200 

(Gey] 



1600 



2000 



1600 



2000 



FIG. 7. The SUSY-QCD correction Aa/a^' and a/asM versus Ms in Case B. The solid curves are for the LHC and the 
dashed for the Tevatron. For Aa/cr", each solid curve overlaps with the corresponding dashed one due to the tiny difference. 



(3) Case C: Only the gluino mass gets much larger than other SUSY parameters (collectively denoted as A/5) ^: 

m~g > Mqjj^^ ~ At,b fJ-r^ Ma-- Ms. (4.5) 
In this case the SUSY-QCD correction behaves as 



Aa 



2as 
37r 



2Ms 



(tan/3 + cota)(l - log — |) 



M? Ms cot a ml 



Ml 



3Mg 



Ml 



Ms tan P m| cos /3 sin {a + (i) ml tan^ fi cot < 



Ml 



sma 



(4.6) 



The main character of this case is that as gluino mass gets large, the correction drops very slowly like log , 

which was also observed in Refs. [11,12]. Again, like other cases, the size of the correction is enhanced by large tan/3. 
In Fig. 8 we show the dependence of the SUSY-QCD correction and a/asM on the gluino mass. Note that in this 



^In this case the Higgs mass bound requires Ms to be much larger than electroweak scale. 
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case we found that tan/3 = 8 cannot satisfy the experimental bounds in Eq.(3.10) for Ms = 600 GeV. In Fig. 8 the 
correction size is significantly smaller than those in Case A and B. The reason is that here a large Ma is chosen so 
that tan/3 + cos a is suppressed (see footnote 3). 

Let us explain the origin of the slowness of the decoupling in this case. Such slowness of the decoupling arises from 



the first term in Eq.(4.6), i.e., -|p(tan /? + cot a) log (note that the AIs in the factor is ^ and the one in the 

logarithm is squark mass Mq). As Mg gets much larger than Mg and /i, decreases but log ^ increases. For the 

600 GeV, the factor is 



example shown in Fig. 8, i.e., Mg is changing from 1 TeV to 5 TeV with fixed Mg 



Mi 



decreased by 1/5 but the factor log ^ is increased by 4.16. Thus the slowness of the decoupling as gluino gets heavy 
is caused by the enlarged mass splitting between gluino and squark. Of course, since tan (5 + cot a is proportional to 
tan /3 cos 2/3 (see footnote 3), the contribution of the first term in Eq.(4.6) will be decoupled rapidly if Ma gets 
large. 
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FIG. 8. The SUSY-QCD correction Acr/o-" and a jasM versus the gluino mass in Case C. The solid curves are for the LHC 
and the dashed for the Tevatron. For Aa/o-", each solid curve overlaps with the corresponding dashed one due to the tiny 
difference. 

(4) Case D: One of the sbottoms and Af ^ become heavy while other mass parameters (denoted as M) are fixed. 
We choose 



or equally 



M^ ~ Mjj ~ At.b Mq ^ mg ^ fi ^ Ma ~ M > Mew 



^ At^b > rn^^ ~ m^^ ^ rUg ^ fj, ^ M :s> Mew, 



(4.7) 
(4.8) 
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where m^^ ^ and m^^ ^ are the masses of bottom-squarks and top-squarks, respectively. 
In this case the SUSY-QCD correction behaves as 



Aa 
'7^ 



2a, 



2M2 ^ / M2 

— — (tan p + cot a) 1 + log — g" 



b2 



f>2 / 



m| cos/3sin(a + /3) 2 2 w^^b^ 



sma 



tan/3) 



(4.9) 



The main feature of this case is the correction decouples like log -^r- and this decoupling behavior is slowed down 

^2 ''2 

by large tan/3. Fig. 9 shows the dependence of the correction and the cross section on (~ ™&2'^' Clomparing with 
the results in case B, we see that the correction size decrease more slowly as gets heavy. 
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FIG. 9. The SUSY-QCD correction Aa/a° and a/asAi versus Alfj in Case D for the LHC (solid curves) and the Tevatron 
(dashed curves). For Aa/cr", each solid curve overlaps with the corresponding dashed one due to the tiny difference. 



From the above analyses we see that when Ma is fixed and all other SUSY mass parameters get large, the SUSY- 
QCD lefts over some remnant effects in the Higgs production process pp ^ bh + X at the hadron colliders. Note that 
for the remnant effects to be left over, fi and gluino mass must be comparable with or larger than the masses of the 
sbottoms. The fundamental reason for such a behavior is that the couplings like hbibj are proportional to SUSY mass 
parameters. 

We conclude this section by making a few remarks. Firstly, in our analysis we assumed mg,fi > and, as a result, 
all four cases have negative values of the correction. In the anomaly- mediated SUSY breaking scenario [30] , a negative 
nig is predicted and in this case, the sign of the correction may be reversed. Secondly, it should be noted that in 
the calculations of such Higgs processes it is necessary to use the loop-corrected relations of the Higgs masses and 
the mixing angle since such relations can significantly affect the results. Thirdly, in case of a light Ma, although the 
SUSY-QCD corrections tend to reduce the cross section severely, the MSSM cross section can still be several times 
larger than the SM prediction. If the cross section of this process is measured in the future and found to be several 
times larger than the SM prediction, a light Ma is favored. 
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V. CONCLUSION 

In this work, we studied the one-loop SUSY-QCD quantum effects in the Higgs production pp [pp) hh + X aX the 
Tcvatron and the LHC in the framework of the MSSM. Wc found that for a hght Ma and large tan /?, the corrections 
can be quite sizable and cannot be neglected. We performed a detailed analysis on the behaviors of the corrections 
in the limits of heavy SUSY masses and found that when Ma is fixed and all other SUSY mass parameters get large, 
the SUSY-QCD lefts over some remnant effects in the Ifiggs production process pp ^ bh + X. Such remnant effects 
can be so sizable for a light Ma that they might be observable in the future experiment. The exploration of such 
remnant effects is important for probing SUSY, especially in case that the sparticles are too heavy (above TeV) to be 
directly discovered in future experiments. 
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APPENDIX A: EXPRESSIONS OF FORM FACTORS 



Before presenting the explicit form of CjS, we define the following abbreviations: 

s = {pi + kf, i={k-p2f, (Al) 

ai,2 = '^(^^^^'> T cos^b), 61,2 = -^(cos^fe ± sin^fc), (A2) 

Af = {ai-bif, Af={ai + bi)^, Aj = aj-bj, (1 = 1,2), (A3) 

Afj = aiaj + bib J + aibj + bjaj, Afj = ajaj + bjbj - ajbj - bjaj, (A4) 

Bjj = aiaj - bib J - ajbj + bjaj, Bfj = ajaj - bjbj + ajbj - bjaj, (A5) 

cci = ^^^^ (7:- I sin^ 9w]sm{a + (3), ( A6) 

cosuw ^ / 

CC2 = "^f ^ sin^ sin(a + /3), 003 = - — — — -( Ab sin a + iJ, cos a), (A7) 

cos Ow o 2mw cos p 

Qii = ccl cof? 6b + cc2 sin^ 9b + 2cc3 sin 9b cos 9b, (A8) 

Q12 = (cc2 - ccl) sin 9b cos 9b + cc3(cos^ 9b - sin^ 9b), (A9) 

Q21 = (cc2 — ccl) sin 9b cos 9b + cc3(cos^ 9b — sin^ 9b), (AlO) 

Q22 = ccl sin^ 9b + cc2 cos^ 9b — 2cc3 sin 9b cos 6b, (All) 

-B/ = -Bi(p,m§,mj J |p2=„2, (A12) 

Bt ^ = Bi{p,mg,mi^)\p2^s, Bj^ = Bi{p,mg,mi^)\p2^i, (A13) 

Cij' = Cij{Pi,k,mi^,mg,mg), Cfj^ = Cij{-pi,-k,mg,mi^,mi^), (A14) 

Cfj ^ = Cij {-p2, k, , m-g, m-g), Cfj ^ = Cij {-p2, k, m-g, m^^ , m^J, (A15) 

C?. " = dj {-p2 , -Ph , m-g , ,mij), C{^" = Cij (-pi ,Ph,mg,mi^,mi^), ( A16) 

Dfj" = Dij{-pi,-k,ph,mg,mij,mi^,mi^), = Dij{-pi,ph,P2,mg,mi^,mi^,mg), (A17) 

Dij" = Dij{-p2,k,-ph,m~g,mbj,mbj,mij), (A18) 

where Bi, Cij and Dij are loop functions defined in [31]. 

2 

After is factored out, the renormalization constant of b-quark can be expressed as 

2 2 

5Zl = CfJ2 ^i^i' = 5^ AfB{, (A19) 
7=1 7=1 
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Smb ^ 



rub ^ mb 

where Cp =4/3. The contributions of the self-energy diagrams of 6-quark propagator can be written as 

2 2 



1=1 1=1 

2 2 

^% = CfY1 AfBt ' - 6Zn, = ^ AfBl ' - 5Zr. 



1=1 1=1 
Ci appeared in Eq.(2.5) are given by 

2 



2 o 1 

= E{- o(^'^?2' + ^^23' + - 1/2 - m?Co« 0^7 + ^^^^2/^7/^ 



1=1 



|(iCr/ + iC^V + - 1/2 - m?Co= /t + ^C|4'^f /*} 



-it/f /lcT7 .Ifr^ . <^"^&N/1 l^ 



-JZi/s - 5Znlt - El/s - nyt - {-5Zl + -5Zr + —^){- + 

I z mb s t 



Tfv}) Sill ^ ^ -j^ 
2 o 1 

^4 = 9(^^?2' + ^^23' + 2CI4' - 1/2 - m?Co« 0^?/^ + ^C72V^?/^ 



7=1 



|(tCiV + tC^V + 2C2V - 1/2 - m?CS ^)^^^/t + Ict.'A'l/i} 



-^t il As: 17 .If'^ , <^"T-b\/l , 1- 



-^Zh/s - 5ZL/t - Sfj/s - Ei/t - {-5Zl + -5Zr + -^){- + -) 

I I rrib s t 

2 

2mw_cos^ 'sr^ A jj 4 . u ^ . 3 ^ ^ 

^n^"sina~ ^^■^'■3^ ° "' 3^ ° ^iJ'*+2 S ° 

2 O 1 

^5 = E{-9(-4C24' + 1 + 2m?Co" - -(sCiV + sC^V + 2C'^M^/s} 

1=1 ^ 

2 

+26Zl/s + 2I]l/g + (JZ^ + 6Zn + 2^)/g + V Q7j{-^mgCo^ "S.Vs 

TOft mfeSina j'^^^ 3 



= E{-|(-4C2«4' + 1 + 2m?Co« ')Af/s - ^(aCf/ + ^C^V + 2C72V)^f /^j 



2 

2 

+2^ZK/g + 2S|,/5 + i5Zr. + JZ^ + 2^)/g + ^"^'^ ^ E Qij{~msC^o " Bfj/s 

TUb mbsma fj^^ ^ 

2 

2 



<^9 = E{-|(-4C2V + 1 + 2m?CS ^)Af /£ - ^(iCf/ + iC^,' + 2Ci/)Af/i} 



1=1 



+2SZn/i + 2Eyi + (JZ^ + dZn + 2^)/f + ^^'^ ""^ ^ V Q77{-^m,C/ "B,Vi 
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Cio = E{-^(-4C|/ + 1 + 2m?CS ')A^li- i(£(7f/ + tCt^' + 2Ci/)Af /{} 



(A27) 



7=1 




(A28) 



Since we have neglect the b-quark mass throughout this paper, Ci, 2, 7,8, 11,12 are irrelevant to our result and we do 
not present their explicit forms here. 
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